Third Neighbor Correlators of Spin-1/2 Heisenberg Antiferromagnet 



Kazumitsu Sakai, 1 Masahiro Shiroishi, 1 Yoshihiro Nishiyama, 2 and Minoru Takahashi 1 

'institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 277-8581, Japan 
2 Department of Physics, Faculty of Science, Okayama University, Okayama 700-8530, Japan 

(Dated: February 27, 2003) 

We exactly evaluate the third neighbor correlator (SjSj +3 ) and all the possible non-zero cor- 
relators (Sf S? +1 S] +2 S$ +3 ) of the spin-1/2 Heisenberg XXX antiferromagnet in the ground state 
without magnetic field. All the correlators are expressed in terms of certain combinations of loga- 
rithm In 2, the Riemann zeta function C(3), C(5) with rational coefficients. The results accurately 
coincide with the numerical ones obtained by the density-matrix renormalization group method and 
the numerical diagonalization. 
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Since Bethe's pioneering work of the spin-1/2 Heisen- 
berg XXX magnet [1], 
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exact evaluation of the correlation functions has been a 
long standing problem in mathematical physics. Espe- 
cially significant are the spin-spin correlators (SjSj +k ) 
(or equivalently (S+Sj +k )/2- Sf = Sf ± iSj), for which 
only the first and second neighbor (k = 1,2) have been 
calculated so far: 



(SfSf +1 ) = — - - In 2 ~ -0.14771572685 
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| . ) p . • - 0.06067976996, (3) 

where £(s) is the Riemann zeta function and (• • ■) denotes 
the ground state expectation value of the antiferromag- 
netic model (J > 0). Here we have taken the thermo- 
dynamic limit L —>■ co. In this letter, we would like to 
report our new results regarding the third neighbor cor- 
relators. Our main result is 

{SjSj +3 } = -^(s^s^^) 
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- yj- C(5) + y C(5) In 2 ~ -0.05024862726. (4) 

In addition, we obtain the third neighbor one-particle 
Green function (cjCj+s)i 

(c]c J - +3 ) £ = i - 2 In 2 + ^((3) - ^C(3)ln2 
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0.08228771669, (5) 



for the isotropic spinless fermion model corresponding to 
(1) by the Jordan- Wigner transformation: 
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Here (...){ denotes the expectation value in the half- 
filled state of the spinless fermion model. Moreover 
we exactly calculate all the possible non-zero correlators 
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The result (2) comes from the ground state energy of 
(1) derived by Hulthen in 1938 [2]. The result (3) was 
obtained by one of the authors in 1977 [3, 4] via the 
strong coupling expansion for the ground state energy of 
the half-filled Hubbard model. This result is also repro- 
duced in the framework of the asymptotic Bethe ansatz 
for an integrable spin chain with variable range exchange 
[5]. However, probably due to the complexity of the wave 
function for these models, no one has succeeded in gen- 
eralizing the method to obtain the higher neighbor cor- 
relators. 

On the other hand, utilizing the representation theory 
of the quantum affine algebra U q {sl2) and the associated 
vertex operators, in 1992, Jimbo et al. derived a univer- 
sal multiple integral representation of arbitrary correla- 
tors for the massive XX Z antiferromagnet [6, 7]. Their 
result has been extended to the XXX [8, 9], the massless 
XXZ [10, 11] and the XY Z [12] antiferromagnets. How- 
ever the explicit evaluation, even for the second neighbor 
correlator (3), was not achieved for a long time. 

In this respect, it is quite remarkable that Boos and 
Korepin recently devised a general method to evalu- 
ate the multiple integral representation especially in the 
study of the Emptiness Formation Probability (EFP) for 
the XXX antiferromagnet [13, 14]. The EFP, P(n) de- 
scribes the probability of finding a ferromagnetic string 
of length n in the antiferromagnetic ground state [9] . Ex- 
plicitly it reads 



P(n) = (n 
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By reducing the integrand of the multiple integral repre- 
sentation to certain canonical form, the EFP for n = 3, 4 
[13, 14] and n = 5 [15] were evaluated by Boos et al. 
(see also recent progress for n = 6 [16]). Note that 
P(2) and P(3) in (7) are related to the first and sec- 
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ond neighbor correlators as P(2) = 1/4 + and 
P(3) = l/8+(5|5| +1 ) + (5^ +2 )/2. 

Here we quote the explicit form of P(4) obtained in 
[13, 14], which is closely related to the third neighbor 
correlator (S?SJ +3 ). 
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P{4) =16 + 4 W+i> + 2 ( ^ +2> + 4< S ^> 
+ (SjSj +1 Sj +2 Sj +3 ) 
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Note that on the antiferromagnetic ground state with- 
out magnetic field, all the correlators with an odd 
number of S z vanishes. Substituting (2) and (3) into 
(8), one finds the relation between the third neigh- 
bor correlator {S z S z +3 ) and the four point correlator 
(S?S? +1 S!! +2 Sf +3 ). However the exact value of (S?S? +3 ) 
itself can not be determined solely from P(4). 

To determine (SjSj +3 ), we consider the following aux- 
iliary correlator: 
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+ (SjSj +1 Sj +2 Sj +3 ). 



(9) 



Here and hereafter P^^et ( a l so written as P^ for sim- 
plicity) denotes a correlator of the form 
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where Sj, ij = {+, — } and E^ 3 is the 2x2 elementary 
matrix ( J E ±± = ±S Z + 1/2, E~+ = S+, E+- = S~) act- 
ing on jth site. In this notation P(4) = Pffff. Note 
that because the Hamiltonian (1) has the symmetry un- 
der the group SU(2), P£ possesses a property like 

if = pi = pi! = pzi (ii) 

As in the case of P(4), the correlators (10) enjoy the 
multiple integral representation [6, 7, 8, 10]: 
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^t/(A 1 ,...,A 4 )T(A 1 ,...,A 4 ), (12) 



where the integration contour C is taken to be a line 
[— oo — ia, oo — ia\ (0 < a < 1). For convenience, we 
choose a = 1/2. The integrand U(X\, . . . , A4) is given by 



U(Xl, . . . , A 4 ) = 7T 
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while T(Ai,...,A4) depends on the selection of s and 
s. Here and hereafter we use the notation Xjk = Xj — 



Xk to save space. In particular for the correlator (9), 
T(Ai, . . . , A 4 ) is given by 

m x ^ A 1 (A 1 +i) 2 A|A 3 (A 3 + ») 2 A| 

1 (Ai, . . . , A 4 ) = — — — - — - — - — — -. (14) 

{A21 — «JA3iA 4 iA32A 4 2(,A 4 3 — I) 

To calculate the multiple integral (12), we follow the 
method by Boos and Korepin [14]. Roughly, their 
method is described as follows. First Taking carefully 
into account the property of U{X\, . . . , A4), we modify 
the integrand T(Ai, . . . , A 4 ) such that the integral gives 
the same result as the original one ("weak equivalence"). 
In this way it is likely that the integrand T(Ai, . . . , A 4 ) 
can always be reduced to the following form (we call it 
"canonical form"): 

T c =P (A 2 , A 3 , A 4 ) + Pl(A \' A3 ' A4) + ^ft^, (15) 

A21 A2lA 43 

where Pq, Pi and P2 arc certain polynomials. Once one 
derives the canonical form, one can perform the multiple 
integral (12) by using the Cauchy theorem [14]. The re- 
sult is written as combinations of the logarithm In 2, the 
Riemann zeta function ((3) and ((5) and rational num- 
bers. Consequently the main part of the calculation for 
the multiple integral (12) reduces to finding the canonical 
form (15). 

Now we consider the case (14) and show that the Boos 
Korepin method is also applicable to our case. Let us 
introduce the following diagram. 
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where j > k. First we expand (14) through partial frac- 
tions. The result consists of 24 terms. Taking into ac- 
count the antisymmetry of the function (13) under trans- 
position of any two variables Xj and A& and the symmetry 
of T (14) under A 4 <-» A 3 , A 2 <-> A 4 , we can reduce the 
24 partial fractions to eight ones. Diagrammatically, its 
denominators are written as 
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From a symmetry of the denominator, the second term 
can further be simplified as 
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where g denotes the numerator of (14). Next using the 
Cauchy theorem, we shift variables Xj — ► \ } ■ ± i such that 
the denominators do not contain i. For instance, we have 




where 

g {1) = -X 1 (X 1 +i)XlX 3 {X 3 + i) 2 Xl 
= - Xl (X 1 +i)X 4 2 X 3 (X 3 +i) 2 Xl 

Finally again using the antisymmetric property of (13), 
we eliminate the symmetric part with respect to Xj <-> A^ . 
In this way we have obtain the canonical form as 
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+ Ai(-4iA 4 + 7Xl - 32iX 3 Xj - 10iX 3 4 - 12iXjXl) 
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Pi =^A 4 - iX\ + ^A 3 A 2 + ^X\ + 2iX 3 X\ + 4A 2 A 4 1 



A?(4A 4 - m\\ - 28X 3 X 2 4 - 10Xl 



- 28tX 3 Xi - S2XiXi), 
P> - — + ;; ; A, + ^AxA 3 - ^A 2 + AiX^j + 6A 2 A 2 . 
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Subsequently, applying the method developed in [14] , we 
calculate the multiple integral by substituting the above 
canonical form into (12). Explicitly, 

pizii = Jo + .h + -h 

J ° = W * = -§ + ^«3) + fjc<5), 
^-4c(3) + ic(3)ln2 + A C (3) 2 -|c(5) 



-C(5)ln2, 
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where Jk denotes the result of the integration regard- 
ing the term P&. We remark that the canonical form is 
not unique due to the non-uniqueness of partial fraction 
expansions. Accordingly, the explicit value of each Jk 
depends on the choice of the canonical form. The final 
result Jo + J\ + J2, however, is always unique as a matter 
of course. 

Combining the result (16) with (8) and (9), we obtain 
the third neighbor correlator (SjSj +3 ) (4) and at the 
same time the correlator (S? Sj +1 Sj +2 Sj +3 ) as 
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W+iWW = go - 3 ln 2 + 30 C(3) 3 C(3) ln 2 
^C(3) 2 -Sc(5) + Sc(5)ln2. (17) 



Now let us consider the four-point correlators of the 
form (SfSf +1 S] +2 S' j+3 ) ({a, £7,*} € {x,y,z,0}; S° = 
1). Because the correlators with an odd number of 
g{<*}\0 vanjgi^ the possible non-zero correlators are re- 
stricted to the following three types: (Sf Sf +1 S? +2 S? +3 ) , 

(S«sP +1 S« +2 SP +3 ) and (SfSf +1 Sf +2 Sf +3 ). Further due 
to the isotropy of the Hamiltonian (1), one can find that 
the independent correlators are written as the follow- 
ing seven ones: (Sf S* +1 S] +2 Sf +3 ), (S* Sf +1 S* +2 Sf +3 ) , 
(SjSj +1 Sj +2 Sj +3 ) and already obtained ones in (2)-(4) 
and (17). Then we shall calculate the remaining three 
correlators here. For convenience we use the operator 
S ± (^ S x ± iS v ) instead of S x . First we consider the 
correlator (Sf S* +1 S* +2 Sr +3 )(= 2<SJS| +1 S| +2 S^ 3 »- 
From (10) and the property (11), this correlator is ex- 
pressed as (S+ S* +1 S* +2 Sr +3 ) = (P++++ - P++r+ ) /2. 
Using the relation (S+ Sj~ +3 ) = 2(P++++ + P+|r+ ), we 
obtain 

(SfSj +1 Sj +2 S j+3 ) = P++++ - -^{SjSj +3 ). (18) 

Here we have used the relation (Sj~Sj~ +3 ) = 2(SjSj +3 ). 
Similarly the other correlators are given by 

/ QZ Q + QZ Q— \ _ p + + + - ^ / QZ QZ \ 

\ a 3 a j+l°j+2 a j+3/ ~- f +-++ ^ j j+2 '' 
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Therefore our goal is to evaluate the auxiliary correla- 
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P+t+ + and P++ 't + ■ They are given if we 



tors P_ 

replace the integrand T(Ai, . . . , A 4 ) by 



T( i) = (Ax + z) 3 A 2 (A 2 + Q 2 (A 3 + z)A§A^-'(A 4 + i) 1 ' 1 
(A21 - i){X 3 i - i)(X 32 - i)A 4 iA 42 A 43 

in the multiple integral representation. Here the correla- 
tor P++++, P+++- and P+++ + correspond to I = 1, 2 
and 3, respectively. Using the procedure similar to the 
case of P+~+~, one obtains the explicit values of the 
above auxiliary correlators. As a result, combining the 
identity (18) and (19) with (2)— (4), we arrive at 
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--C(3) In2--C(3) 2 - -C(5) + -C(5) In 2, (20) 

(^sr« +1 sr J - +a s» +fl ) = I L-Iin2 + ^c(3) 

5 3 65 5 

--C(3)ln2--C(3) 2 --C(5) + -C(5)ln2, (21) 



(SfS j+1 Sj +2 Sj +3 ) 
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^ + ^l n2 -_C(3) 
120 6 120 v ' 
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+ -C(3) ln 2 + -C(3) 2 + -C(5) - -C(5) In 2. (22) 

We mention a few remarks of our results, (i) All the 
above correlators are written as the logarithm ln2, the 
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TABLE I: Estimates of the correlators by the exact evalua- 
tions, DMRG and the extrapolations from the numerical di- 
agonalization for the system size L — 24, 28, 32. 



Correlators 



Exact 



DMRG 



Extrap. 



I Q + QZ QZ q— \ a 

Wj ^+1^+2^+3/ 

/ QZ QZ QZ QZ \ 

(Sj Sj +l Sj +2 Sj +3 ) 



-0.0502486 
0.0205719 
0.0307153 

-0.0141607 
0.0550194 



-0.0502426 
0.0205681 
0.0307105 

-0.0141579 
0.0550108 



-0.0502475 
0.0205716 
0.0307154 

-0.0141606 
0.0550198 



Riemann zeta function £(3) and ((5). This agrees with 
the general conjecture by Boos and Korepin: arbitrary 
correlators of the XXX antiferromagnet are described as 
certain combinations of logarithm In 2, the Riemann zeta 
function with odd arguments and rational coefficients. 
Especially intriguing is the existence of the non-linear 
terms such as ((3) 2 , ((3) In 2 and ((5) In 2. (ii) The corre- 
lator (20) is interpreted as the third neighbor one-particle 
Green functions (c]cj +3 )f /4 via the Jordan- Wigner trans- 
formation (6). Obviously, the first neighbor one-particle 
Green function is expressed as 



<4 c i+i)f 



In 2 



-0.295431453707, 



which coincides with (Sj~ Sjj_i ). Due to the characteristic 
(cjcfe)f ~ A jk cos(ir(k - j + l)/2), one 
= 0. Therefore the quantity (5) is the 



J 3 ^3+1/ 

7r/2-oscillation: 

finds (c]cj +2 fc}f 
first non-trivial exact result of the correlators containing 
the fermionic nature, (iii) The difference between (20) 
and (21) gives the nearest chiral correlator 

((SjxS j+1 ) ■ (S j+2 x S j+3 )) = 3((S+S* +1 S- +2 S* +3 ) 



(SfSj +1 Sj +2 S j+3 )) 



which exactly agrees with the one derived from the 
ground state energy of an intcgrable two-chain model 
with four-body interactions [17]. 

To confirm the validity of our formulae, we performed 
numerical calculations by using the density-matrix renor- 
malization group (DMRG) [18, 19] and numerical diag- 
onalization. As for the DMRG, wc followed standard 
algorithm [20]. We have repeated renormalization 500- 
times. At each renormalization, we kept, at most, 200 
relevant states for a (new) block. The numerical diago- 
nalization was performed for the system size L = 24, 28 
and 32. We extrapolate the data from a fitting function 
a + ai/L 2 + a 2 /L 4 . All our analytical results coincide 
quite accurately with both numerical ones (TABLE I) . 

In closing we would like to comment on generalizations 
of the present results. The extension to the calculation of 
higher neighbor correlators (S* Sj +k )k>4 is of great inter- 
est. The fourth neighbor one (S , |S , | +4 ), for example, will 



be calculated by combination of the EFP, P(5) and two 
independent auxiliary correlators, which can in principle 
be evaluated. In fact P(5) has been already obtained in 
[15]. The computation, however, will be much more com- 
plicated. Alternatively, extending the present result to 
the inhomogeneous case as in [16] and taking into account 
the property of the quantum Knizhnik-Zamolodchikow 
equation, we may derive higher neighbor correlators. Us- 
ing this, eventually we hope to extract the long-distance 
asymptotics (<SJ<S , j : +fc }fc»i, which is a crucial problem in 
conformal field theory [21, 22]. 
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